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^^ ' Abstract 

r-{ , The approximation for the NJL gap equation that was developed in 

our previous paper allows us to investigate vacuum inhomogeneities in 
the mean field approach. The simplest case of a domain wall is studied 
thoroughly. The Jackiw-Rebbi problem which arises in treatment of the 
fermions' interaction with the domain wall is solved explicitly and the 
entire spectrum of bound states is described. It is estabilished that under 
our assumptions no bound states exist at the domain wall, however higher- 
order corrections may make them appear. Thus existence of the quark 

Mh, bound states at domain walls may be a good test for subtle features of 

^ . QCD. 

1 Introduction 

>, 

0^ , Deriving quark confinement from QCD is one of the most intriguing problems 

in modern quantum field theory. Nonlinearity of theory's Yang-Mills equations 
is a huge obstacle. It's well known that at the low energies spontaneous gauge 
symmetry breaking occurs. However the gluon propagator structure is unknown 
which makes it difficult to build any plausible approximations. 

Nambu-Jona-Lasinio model first proposed long before QCD was developed[l] 
is believed to be one of such approximations. In one of our previous works [2] 
we've proposed a derivation of it from QCD in the mean field approach to 
the latter developed by Kondo|3] It's notable that in our approach NJL model 
is obtained naturally together with a cutoff parameter which is equal to the 
effective gluon mass. 



A = niA (1) 

Another nontrivial feature of our approach is that the interaction constant 
of the effective NJL lagrangian depends on the cutoff parameter too: 
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p is a free parameter that arises in the mean field approximation to QCD 
and N is the number of colors which we are keeping for the sake of possible 
generalisations. 

Mean field approach is also a paradigm for dealing with NJL model itself. 
It was first applied by Nambu and Jona-Lasinio in their seminal work where 
they've found the model's condensate assuming homogeneity i.e. that the mean 
field is constant over the spacetime. Finding inhomogeneous configurations 
is tricky since one needs to solve the nonlinear integral gap equation. Little 
progress has been made on this path so far with the works of M. Thies et al. [3] 
being a notable exception. 

We've also followed this path in our last work[5' . There, we've developed a 
differential approximation for the gap equation based on the study of divergen- 
cies of the mean field action. The equation turned out to be Landau- Ginzburg- 
like containing derivatives up to second and nonlinearities up to third order. 
But here the superconductivity analogy ends since the equation is much more 
complicated than plain Ginzburg-Landau. Even the mean field itself is much 
more complex being a matrix variable: 

^^S. + in'' + v + j^iu. (3) 

However we've found that this equation has a simple scalar sector. In other 
words it allows to put 17 = ^; ^ G R. The equation then gets much more simple 
and if the NJL model fermions are initially massless it gets reduced to the 



Dax)~Kax)+2eix)=0, (4) 

where K is a constant that depends on the model's parameters and is quadratic 
in the cutoff K ^ A^. 

This equation has the famous one-dimensional kink solution 

^ (a;'') = fi tanh [Xz] ; z = x^ (5) 

which can be interpreted as a domain wall between two vacua. It's worth noting 
that the exact gap equation of NJL model has two solutions with opposite signs 
in the homogeneous case so the domain wall configuration is not an artifact of 
our approximation. 

The mean field approach to NJL model suggests that after the mean field 
equations have been solved we should substitute the solution in the Dirac equa- 
tion to find the cortresponding fermionic states: 

id - hixf")] ij ^ (6) 

If we substitute the kink solution (IS|) for ft we get the well-known Jackiw-Rebbi 
problem [6j. Jackiw and Rebbi investigated the interaction of Dirac fermions 
with a 4>'^ scalar field kink and found the zero-energy localized fermionic state. 
Our goal here is to pick up the task of solving ([6]) completely filling the gap 
in their work. It's also worth noting that while the original Jackiw-Rebbi case 
of fermions interacting with a real scalar field is somewhat artificial our NJL 
background gives this problem an interesting physical interpretation. 



2 Jackiw-Rebbi problem in the context of NJL 
model 



So our task is to solve the equation 

id — /i tanh (Az) 
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But let us first find the meaning of the kink parameters in the context of NJL 
model. 

In [5] we've found that the exact form of the 0^ equation Q is 



where 



and 



nC(x) - [2 - G(7V, p)] Z(A, M)ax) + 2^3 (^) = q 

32^2 (p2 _ 1^ 



(8) 



G{N,p) 



Z{A, M) = 



N 
A2 



In 






1 



A is the QCD cutoff parameter and p > 1 is a free parameter that arises 
in the mean field approach to QCD. The meaning of p is unclear by now and 
revealing it can be a task for further investigation. 

M is found from the homogeneous gap equation 
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Having the homogeneous gap equation, we can notice that the (j)^ equation 
now becomes 



n^(x) - 2M^i{x) + 2£,^{x) = 



(10) 



An interesting case here is the large M case M ^ A. After expanding the 
logarithm and neglecting the higher order terms the gap equation turns then 
into 



M^ = 



NK^ 
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and we can conclude that the large M approximation is valid when p — >■ 1 . 
A kink solution for the equation ([TU| is 
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M tanh (Mx'^ 
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3 Solution 

Let us now solve the Dirac equation with 0^ kink potential. It's enough to study 
the special case when the equation becomes 



[i-f°dr + i^^d^ - p tanh (z)] -0(r, z) = 



(13) 



Eigenstates that depend on other spatial coordinates can be constructed from its 
eigenstates by the means of Lorentz transformations. The "wideness" parameter 
of a kink can be absorbed into its "height" by scaUng the spacetime variables 
appropriately: 






(14) 



The next step is putting 
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By choosing the appropriate Dirac matrix representation we reduce the equation 
to the following system: 



dx 



iaz h [-E — /itanh(z)] = 

dz 



d<j) 



iaz— — \-[E + fita.iih{z)]x — 
dz 



Then we make an ansatz 



0(z) = 0(z)|t) 
X{z)^x{z)n) 



(15) 



(16) 



where cr^ |t) = |t) (we might have chosen the other eigenvector of the Pauli 
matrix which wouldn't change the picture much). 

The resulting system of equations for (f>, x can be written in the following 
matrix form: 

[iE(Tx — /i tanh(z)crj;j 



dz 
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(17) 



Now let us define 
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These vectors satisfy the following identities: 



We can now make the following ansatz: 

== p(z)cosh"^(z) 
X 



+ q{z) cosh^iz)\-] 
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(20) 



which leads us to 
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(21) 

Then we exclude q{z) from this system of equations and obtain the following: 

cosh^ ''(z) dp 



q{^) 



E dz 



4| - 2^ tanh(0) ^ + E^p{z) = (22) 

dz^ dz 

If i? = p and q decouple from each other and we can just put p — 1,9 = 0; the 
other linearly independent solution diverges at the infinities and we don't take 
it into account. 

Finally, we put Q = sinh(z) and the last equation turns into 

(1 + C') |l + (1 - 2m) c| + EMO = (23) 

This equation falls into the hypergeometric class, furthermore, it can be eas- 
ily proven that there's a series of polynomial solutions that are orthogonal at 
(— oo, oo) with the measure 

W{0 = {1 + CT'^'^. (24) 

Their spectrum is 

E^„^n{2n-n) (25) 

It can be obtained by studying the asymptotics of the equation and taking into 
account that the leading order should disappear. These solutions with \n\ < ^ 
correspond to the bound states. The Rodrigues formula for these polynomials 

is 

p^{o = {i+eY^'^^^{i+e)-'-^^. (26) 

One can note that the polynomials are very similar to the well-known Gegen- 
bauer polynomials. See the appendix for more details. 



4 Discussion 

Now, we have established that the equation ([T3| has a series of bound-state 
solutions with the spectrum (^5]) . However the E = Q solution that does not 
have its opposite-energy counterpart should be considered the fermionic vacuum 
and the symmetry between opposite energy solutions corresponds to fermion- 
antifermion dualism. Therefore, we can conclude that the actual bound states 
exist only if /i > 1. 

Let's now recall the transformation formulae ([T^ and the actual parameters 
of the NJL kink p^ . After doing this we can make a very interesting observation 
that under our approximation exactly the "borderline" /i = 1 case takes place. 
Furthermore, this does not depend on the NJL condensate and thus it doesn't 
depend on the cutoff! 



Let us now look once more at the scaling formulae ([T^ . The one of interest 

is 

H 

where H is the height of the kink and j is its wideness. Therefore we can 
conclude that the wider is the kink the deeper the effective potential well be- 
comes. So we may expect that if some subtle effects of QCD effectively widen 
the domain wall then binding of fermions to it can take place. 

So existence of fermionic bound states at a domain wall may be a good test 
for those subtle effects. 



5 Appendix. Orthogonality 

The bound states' wavefunctions should be orthogonal: 

cPx'4>^{x, y, z)il;rn{x, y, z) = N{n)Snm- (27) 

In our case the states are actually "semibound" since the fermion motion 
is restricted along the z-axis along. Thus the normalization constant will be 
infinite. So we actually should prove that this scalar product reduces to 



cPx'4>^^{x,y,z)il)rn{x,y.,z) = / dxdy / dz'4>^{z)'4>m{z) 



(28) 



and the z-integral is actually finite and J dzip^ {z)ip„i{z) — CnSnm for the wave- 
functions that don't depend on x, y that were studied above. 

To achieve this we should first reconstruct the complete Dirac spinors for 
the (1261) solutions. Taking into account the transformations from section 3 we 
obtain 



V'™t(^) 



cosh '^{z) 
V2 



(^p„(sinh(z)) + -^ 
i (pn{smh{z)) - ;^ 



dp„(sinh(z)) N j i ^x 
t) 



dz 
dpn{sin.h.{z)) 



(29) 



Thus the normalization condition becomes 



{nt\m^) = I dxdy / dz cosh ^^(z) x 

1 dp„(sinh(z)) (ipm(sinh(z)) 



p„(sinh(z))p,„(sinh(z)) + 



i!m dz dz 

However for the bound states we have, by substituting again sinh(z) = ( 
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Therefore we can conclude that solutions with En = —Em (we can make a 
convention _E_„ — —En,p-n{0 = Pn(C)) are orthogonal. For the rest of the 



+00 



cases scalar product reduces to / d( (^1 + ^^) ^ PniOPmiC) and we can see 



that the measure is exactly ([24|. 

Only the opposite spin case remains now uninvestigated. However orthogo- 
nality of the opposite spin solutions is guaranteed trivially by (t | i) = and 
the spin-down solution are very similar to the spin-up case studied above, that's 
why we aren't going into studying them in detail. 
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